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ABSTRACT 


A bsorption Heat Pumps for Space Applications 
In the first part of this study, the performance of the Absorption Heat Pump (AHP) 
with water-sulfuric acid and water-magnesium chloride as two new refrigerant-absorbent 
fluid pairs was investigated. A model was proposed for the analysis of the new working 
pairs in a heat pump system, subject to different temperature lifts. Computer codes were 
developed to calculate the Coefficient of Performance (COP) of the system with the 
thermodynamic properties of the working fluids obtained from the literature. The study 
shows the potential of water-sulfuric acid as a satisfactory replacement for water-lithium 
bromide in the targeted temperature range. The performance of the AHP using water- 
magnesium chloride as refrigerant-absorbent pair does not compare well with those 
obtained using water-lithium bromide. 

The second part of this study concentrated on the design and testing of a simple 
ElectroHydrodynamic (EHD) Pump. A theoretical design model based on continuum 
electromechanics was analyzed to predict the performance characteristics of the 
electrohydrodynamic (EHD) pump to circulate the fluid in the absorption heat pump. A 
numerical method of solving the governing equations was established to predict the 
velocity profile, pressure - flow rate relationship and efficiency of the pump. The 
predicted operational characteristics of the EHD pump is comparable to that of 
turbomachinery hardware; however, the overall efficiency of the electromagnetic pump is 
much lower. An experimental investigation to verify the numerical results was conducted. 
The pressure - flow rate performance characteristics and overall efficiency of the pump 
obtained experimentally agree well with the theoretical model. 
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Refrigerant -Absorbent Pair Selection for the Absorption Heat Pump 


Refrigerant-Absorbent Pair Selection for the Absorption Heat Pump 


Introduction 

The absorption heat pump (AHP) has come a long way since its invention by 
Ferdinand Carre in 1860. Due to the inherent absence of moving parts, the AHP has the 
advantage of relative simplicity and potential reliability. On the other hand, it requires 
large heat exchangers for acceptable coefficient of performance (COP). It also exhibits 
relatively lower cooling-mode COP as compared to that of Rankine cycle refrigeration 
machinery. The COP of the AHP depends on the following parameters: 

System Configuration The system configuration refers to the various 
components that make up the heat pump. The components of the system are usually 
dependent on the particular application, refrigerant-absorbent pair, and other factors like 
space, weight, etc. The final layout of the heat pump is generally determined by the 
economics of the system. 

Operating Temperature Range . The operating temperature range decides the 
selection and performance of the refrigerant-absorbent fluid pair. An optimal temperature 
range cannot be specified for all applications. The temperature range for an AHP is 
decided by the type of heat source and heat sink. 

Properties of Refrigerant-Absorbent Pair . The thermodynamic properties of 
the refrigerant-absorbent pair are the deciding factor in its selection as a feasible working 
pair. The properties needed to calculate the COP of an AHP have been discussed in detail 
by Alefeld (1987). The desirable properties of the working pairs are described in later 
sections. 

Various system configurations used to improve the COP of the AHP have been 
reported (Serpente, Perez-Bianco, & Seewald, 1992; Garimella, Christensen, & Petty, 
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1992). These modifications have resulted in an improvement in system COP. The 
operating temperatures in almost all practical heat pumps depend on the particular 
application. A temperature range cannot be prescribed for all cases; however, a suitable 
choice of refrigerant-absorbent pair will allow an optimum COP. The most commonly 
used refrigerant-absorbent fluid pairs are ammonia-water (NH 3 -H 2 0) and water-lithium 
bromide (H 2 0-LiBr). The properties and working parameters of these fluid pairs have 
been thoroughly investigated (ASHRAE Handbook, 1985). A few other working pairs 
have been suggested (Bhatt, Krishna Murthy, & Seetharamu, 1992; Kaushik, Gadhi, 
Agarwal, & Kumari, 1988). The basis for the selection of a refrigerant-absorbent pair has 
been described in the literature (Mansoori & Pate, 1979; ASHRAE Handbook, 1985). 
The AHP analysis based on the first and second laws of thermodynamics has also been 
well documented (ASHRAE Handbook, 1985; Stoecker & Jones, 1982). 

In the present effort, an attempt was made to understand the thermodynamics of 
new refrigerant-absorbent pairs. The objectives of this study were twofold. 

1 . To investigate the feasibility of using water-sulfuric acid and water-magnesium 
chloride as new refrigerant-absorbent fluid pairs for the absorption heat pump. 

2. To design an EHD pump for space applications. 

The first part of the study deals with the determination of thermodynamic 
properties of aqueous sulfuric acid and aqueous magnesium chloride solutions. The 
thermodynamic properties of interest in AHP simulation are vapor pressure, solubility, 
relative volatility, latent heat and enthalpy of the fluid pairs. A computer program was 
wntten to calculate the performance of both working fluid pairs in the system. 

The design and development of an EHD pump has been detailed in the second part 
of this study. The results obtained testing the pump are also presented. 
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Characteristics of Refrigerant-Absorbent Fluid Pair 


The efficiency of the absorption heat pump depends on the selection of a proper 
refrigerant-absorbent fluid pair. The operating pressures, heat transfer rates and hence the 
COP of the system are decided by the physical, transport, and thermal properties of the 
refrigerant-absorbent pair. The following section discuss the desirable properties of the 
working pairs. The thermodynamic properties of sulfuric acid solution and aqueous 
magnesium chloride are also described. 

Desirable Properties of Refrigerant-Absorbent Pairs 

Certain chemical and physical properties of the refrigerant-absorbent pair must be 
satisfied before it can be selected for use in the absorption heat pump. The properties of 
interest specified by the American Society of Heating, Refrigerating and Air-Conditioning 
Engineers (ASHRAE) are: 

Absence of Solid Phase . The refrigerant-absorbent pair should not form a solid 
phase over the operating concentration and temperature ranges. The freezing-point 
temperature and crystallization-point temperature of the pair over the operating 
concentration range should be well below the corresponding operating temperature. The 
presence of solid phase depends on the maximum solubility of the refrigerant in the 
absorbent at a given temperature. 

High Volatility Ratio . The volatility ratio is the ratio of the refrigerant volatility 
to that of the absorbent at any given pressure. The volatility ratio directly affects the heat 
transfer rate in the generator. As the volatility ratio increases, the ease of separation 
improves. The volatility of the refrigerant must be higher than that of the absorbent (i.e., 
relative volatility » 1). 
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aren a Affinity . The absorbent must possess strong affinity for the refrigerant 
under normal operating conditions, greater the affinity between the refrigerant-absorbent 
pair the lesser the amount of absorbent to be circulated. 

Mo derate Operating Pressure The operating pressures of the heat pump 
depend on the physical properties of the refrigerant-absorbent pair. High pressures 
necessitate the use of thick-walled equipment while low pressures necessitate the use of 
large volume equipment. Moderate pressures are therefore preferred. 

C hemical Stability The refrigerant-absorbent pair needs to be chemically stable 
under normal operating conditions. Chemical instability can lead to the formation of 
unwanted substances such as solids, gases or corrosives that can damage the equipment. 
Lo w Corrosiveness Almost all the commonly used refrigerant-absorbent pairs 

are corrosive to some extent. Care should be taken to minimize corrosion of mechanical 
parts. 

Safety As far as possible the fluid pairs used in the absorption heat pump should 
be nontoxic and nonflammable. Fluids that are not safe are limited to outdoor use. 

Ir ansport Properties . The refrigerant-absorbent pair should possess such 
properties as low viscosity, moderate density, low surface tension, high thermal diffiisivity 

and high mass diffiisivity. These transport properties have an effect on the overall 
performance of the system. 

High Latent Heat The latent heat of the refrigerant has a large impact on the 

refrigerating effect. The greater the latent heat, the lesser is the circulation rate of both 
refrigerant and absorbent. 

Detailed characteristics of commonly used refrigerant-absorbent pairs have been 
reported at length in the literature (ASHRAE Handbook, 1985). The properties 
considered in the present study are solubility, pressure, enthalpy, relative volatility and 
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freezing point. All the concentrations referred to in this study are with respect to the 
absorbent (H 2 SO 4 or MgCl 2 ). 

Thermodynamic Properties of Water-Sulfuric Acid Fluid Pair 

The thermodynamic properties of interest in the determination of the COP are 
temperature, pressure and concentration. Knowledge of any two of the three implies the 
third. The physical properties of the fluid pair establish the operating pressures of the 
cycle. The properties of water-sulfuric acid solution have been documented in the 
literature (Warrier, Annamalai, Nguyen, & Lin, 1993). The total pressure of the fluid pair 
as a function of temperature and concentration was regressed from data available in the 
literature (Chemical Engineers' Handbook, 1973; CRC Handbook of Chemistry and 
Physics, 1980). The total pressure of the solution as a function of temperature and 
concentration is given in Table 1. The equation in Table 1 is an Antoine-type equation 
with 50 data points utilized for the regression. The percentage error for the regressed 
equation was ±3%. The computed mean square error was 18.5 and the absolute standard 
deviation was 4.3. 

Sulfuric acid is completely soluble in water at all temperatures (infinite solubility). 
In cold water, heat is evolved during mixing (heat of dilution), while in hot water no heat 
is evolved. Since sulfuric acid theoretically has infinite solubility, sulfuric acid solution of 
any concentration may be used. The concentrations used in practice are limited by the 
corrosiveness of the sulfuric acid solution. 

It is essential that the solution does not freeze within the operating temperatures 
and concentrations. The freezing point of water-sulfuric acid solution as a function of the 
solution concentration can be obtained from Chemical Engineers' Handbook (1950). The 
data obtained are plotted in Figure 1 . As the solution freezing points are well below the 
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normal operating temperature limits over most concentrations, the chances for the solution 
freezing and causing equipment shutdown are minimal. 

The ease of separation of the refrigerant from the solution depends upon the 
relative volatility of the refrigerant in the solution at the given temperature and pressure. 
The greater the relative volatility, the easier is the separation. The relative volatility of 
water in an aqueous sulfuric acid solution is determined using the equation outlined in 
Table 2. The relative volatility of the aqueous sulfuric acid solution for various 
temperatures and at the generator pressure is illustrated in Figure 2, where yj and xj are 
obtained from equilibrium vapor pressure data. Figure 3 shows the variation in relative 
volatility with pressure. From Figure 3, it can be seen that as the operating pressure 
increases, the relative volatility increases and the separation becomes easier. This property 
helps provide an initial estimate of the feasible operating pressure range. 

The vapor pressure equilibrium diagram for a mixture is also important during the 
preliminary design stages, since the formation of azeotropes should be avoided. An 
azeotrope is formed when the mole fraction of a substance in a mixture is equal in both the 
liquid and vapor phases (i.e., yj = xj). An azeotrope forms a constant boiling point 
mixture and hence separation becomes impossible. 

The COP of the absorption heat pump depends on the enthalpy of the solution at 
various temperatures and concentrations. The enthalpy for various temperatures and 
concentrations was regressed from enthalpy-concentration diagrams available in the 
literature (Smith & Van Ness, 1979). The enthalpy as a function of temperature and 
concentration is given in Table 3. The effect of pressure on the enthalpy of the solution is 
negligible. The regression was performed using a total of 50 data points. The computed 
mean squared error was 33.14 and the absolute standard deviation was 5.76. The 
percentage error was found to be 4.36%. 
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Thermodynamic Properties o f Water-Magnesium Chloride Fluid Pair 

The thermodynamic properties of aqueous magnesium chloride solution are 
described in this section. Like all salts, magnesium chloride has limited solubility at a 
given temperature. The operating concentration at any given temperature should be less 
than the maximum solubility at that temperature to prevent ciystallization. Crystallization 
of the salt leads to blockage of equipment and shutdown 

of the unit. The maximum solubility of magnesium chloride is given by the correlation 
shown in Table 4. The solubility data was obtained from the literature (Mel'nik & 
Melnikov, 1970). The maximum solubility at various temperatures is plotted in Figure 4. 

As seen from the figure, the solubility limits the operating concentration range at any 
temperature. 

The vapor pressure of water-magnesium chloride solution was obtained from the 
literature (Sako, Hakuta, & Yoshitome, 1985; Patil, Tripathi, Pathak, & Katti, 1991). The 
equations used to calculate vapor pressures of the solution are shown in Tables 5 and 6. 
These correlations in Table 5 and 6 are Antoine-type equations. Available experimental 
data are very limited. No data are available for vapor pressure in the temperature range of 
273K to 303K. Similarly, no vapor pressure data are available for higher molalities. 

Hence, the equations available in the literature were extrapolated to determine the vapor 
pressure. 

The relative volatility of the water (refrigerant) in a water-magnesium chloride 
solution was calculated using the equation described earlier in Table 2. The relative 
volatility is shown in Figure 5 as a function of pressure. The relative volatility of water is 
very high as seen in Figure 5, and hence separation will not be a problem. 

The enthalpy of aqueous magnesium chloride was determined from data available 
in the literature (Snipes, Manly, & Ensor, 1975). Table 7 shows the equation used to fit 
the data points determined experimentally. A total of 43 data points were used for the 
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regression analysis. The mean square error was 0.069 and the absolute standard deviation 
was 0.263. The mean percentage error was calculated to be ±0.04%. Enthalpy at higher 
temperatures were obtained from experimental data (Mayrath & Wood, 1983). The 
regression analysis was done using a total of 46 data points. The mean square error was 
0.002 and the absolute standard deviation was 0.052. The mean percentage error was ± 
0.15%. The correlation equation is shown in Table 8. As seen from Tables 7 and 8, the 
enthalpy correlations do not cover all temperature and molality ranges. The correlation 
equations obtained have been extrapolated to obtain vapor pressure and enthalpy data for 
those ranges for which no experimental data were available. 
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Absorpti on Heat Pump Cycle Analysis 


The properties of the refrigerant-absorbent pairs were discussed in the previous 
sections. This section discusses the analysis of the performance of water-sulfuric acid and 
water-magnesium chloride as new refrigerant-absorbent pairs. The AHP configuration, 

basic governing equations and performance calculation procedure are described in the 
following sections. 

Absorpti on Heat Pump Configuration 

A simple AHP configuration was chosen for this study. The system consists of a 
generator(G), condenser(C), evaporator(E), absorber(A), and a solution-side heat 
exchanger (HE) as shown in Figure 6 (Stoecker & Jones, 1982). The function of a 
compressor in a vapor-compression system is accomplished in a vapor-absorption system 
by the combination of the absorber, solution pump and generator. Heat is supplied at the 
generator to boil off the refrigerant from the strong solution. The generator is a 
separation device where the fluids are separated by distillation. The differences in the 
boiling points between the fluids and the relative volatility of the refrigerant are the factors 
deciding the ease of separation. The temperature and pressure of the refrigerant leaving 
the generator are in equilibrium. The high-pressure refrigerant vapor is condensed at the 
condenser, where heat is released. An expansion device then expands the liquid 
refrigerant to a lower pressure state prior to the evaporator. 

The liquid refrigerant absorbs the heat of vaporization and evaporates in the 
evaporator. The temperature and pressure of the low-pressure vapor leaving the 
evaporator are in equilibrium. This vapor is absorbed at the absorber by the weak solution 
returning from the generator. The heat of dilution is removed in the absorber. The strong 
solution is pumped to the high-pressure generator by a pump. The refrigerant is then 
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regenerated in the generator and the cycle continues. A heat exchanger is provided prior 
to the generator to transfer the heat from the hot, weak solution exiting the generator to 
the cool, strong solution being pumped from the absorber. This heat exchanger has a 
positive impact on the COP of the system as it reduces the supplied heat to the generator 
required to free the absorbate. 

Governing Equations and Performance Calculation 

The following three equations form the basis for the thermodynamic analysis of the 
absorption heat pump cycle. 

Mass Balance . The mass balance equation is given by: 

i>,=o to 

i 

where, m is the mass flow rate. 

Partial Mass Balance (Material Balance) . The conservation of individual mass 
balance is given by: 

i>i£=° (2) 

i 

where, m is the mass flow rate and E, the solution concentration. 

Energy Balance . The conservation of energy is given by: 

ZQ i +E m . | 'i=° (3) 

i i 

where, Q is the heat transfer and m the mass flow rate. The above equations can be 
applied to any state point of the heat pump system. 

The performance calculations are based on the following assumptions: 

1. The load is constant at 1 ton. 

2. Selected evaporator, condenser and absorber temperatures are maintained 

constant. (T a = 303K, T e = 283K, T c = 313K) 

10 




3. Equilibrium conditions exist at the generator and absorber outlets. 

4. The expansion processes are isenthalpic. 

5. The heat exchanger effectiveness (s) is 0.8. 

With a 1-ton load and constant temperatures at the evaporator, condenser and 
absorber, the COP of the system was calculated for various temperature lifts (i.e., 
temperature difference between the hot and cold plates). The condenser/generator 
pressure was the saturation pressure of the refrigerant vapor corresponding to the 
condenser temperature (T c = 303K). Saturation pressures and enthalpies of the 
refrigerant at different state points were obtained from the literature (Raznjevic, 1976). 
The saturation pressure corresponding to the evaporator temperature (T e = 283K) 
dictates the pressure of the evaporator/absorber that forms the low pressure side of the 
system. As the load is maintained constant, the refrigerant mass flow rate is constant. 
The refrigerant mass flow rate corresponding to state point 6 in Figure 6 is given by: 


load _ Q e 

^7 ” ^6 ^7 h 6 


( 4 ) 


where, hg and h 7 are the enthalpies corresponding to saturated water at (T c and P c ) and 
(T e and P e ), respectively, Q e being the heat transfer at the evaporator. 

The following mass flow rate relations are applicable to the system: irq = m 2 , 
m 3 = 1114 , and m 5 = mg = mq. The total mass flow rate of the system is obtained from: 

m, = m 3 + m 5 (5) 

The H 2 SO 4 mass flow rate is calculated from: 

m i 6 = £ (6) 

where, = concentration based on T a and P a and £3 = concentration of solution based 
on Tg and Pg. 

The heat rejected at the condenser and absorber is given by: 
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( 7 ) 


Q c = m 5 (hj -h 6 ) 

Q, = m 7 h 7 +m 4 h 4 -m, h, 


where, m and h are the mass flow rate and enthalpy of corresponding state points. The 
heat supplied at the generator is obtained from the energy balance around the generator: 

Q g = m 3 h 3 +m 5 h 5 -m 3 h 2 (g) 


At the heat exchanger: 


s = 0.8 = 


m . (jvA) 


m, (h, - h;) (9 > 

where, h 3 — enthalpy of solution for £ 3 and T a . The above relation for e is given by Soo 
(1958). 


The specific circulation flow rate is defined as the ratio of the weak solution flow 
rate (mj) to the refrigerant vapor flow rate (m5), as given by: 


F = -Ji±A = -6 

6-6 6-6 


( 10 ) 


since ^5 _ 0 (relative volatility » 1). In Equation 10, F is the specific circulation flow 
rate and £ is the concentration of the aqueous solution. 

Once all state points (1 through 7) are determined, the COP of the system is easily 
calculated. The coefficient of performance (COP) of the system is defined as: 


Heating mode : COP heit 
Cooling mode : COP^, 


= Q. + Qc 

Q, 

Q g 


(11) 

( 12 ) 


where, Q e , Qg, Q c , and Q a are the heats at the evaporator, generator, condenser and 
absorber, respectively. 


The COP of the AHP system depends upon many factors, such as system 

configuration, operating temperature range and fluid properties. The results obtained 
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were compared with those obtained from a heat pump system using water-lithium bromide 
as the refrigerant-absorbent pair. Data for the aqueous LiBr solution were obtained from 
the literature (ASHRAE Handbook, 1985). The above equations were solved for a range 
of temperatures at the generator. Computer codes were developed and the system 
performance analysis was carried out on a microcomputer. The results obtained were 
compared with the vapor equilibrium data for assessment of accuracy. The pump work, 
being negligible, was not included in the calculations for the system performance. 
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Results and Discussion 


The physical and thermodynamic properties of the two new working pairs, a 
description of the heat pump system and the calculation procedures have been covered to 
this point. The results of the absorption heat pump modeling with the newly developed 
working pairs are presented in the following sections. 

Heat Pump Performance 

The coefficient of performance of the heat pump was calculated for both the 
heating and cooling modes. The results obtained for water-sulfuric acid and water- 
magnesium chloride pairs are presented below. 

AHP Performance using the Water-Sulfuric Acid Pair 

The results obtained with water-sulfuric acid as the working fluid pair are shown in 
Table 9. For comparison purposes, results of the water-lithium bromide system are shown 
in Table 10. For the water-sulfuric acid system, the minimum temperature lift is 70K. 
Below this range, the concentration of the strong solution (£ 3 ) is less than that of the 
weak solution (£ 1 ). As a result, the refrigerant mass flow rate (m$) according to Equation 
5 would become negative, which is an impossibility. The material balance around the 
generator thus determines the minimum possible temperature lift . The water-lithium 
bromide heat pump system is also considered for the same temperature lift. 

The cooling-mode COP of the water-sulfuric acid system as shown in Figure 7 , 
increases from 0.788 at T g = 373K, to a maximum of 0.805 at T g = 353K, and then 
decreases rapidly. For the same range, the COP of the water-lithium bromide system 
increases from 0.801 at T g = 373K to 0.817 at T g = 353K, and then decreases. 
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For constant T g , T c , T e , and T a , as the specific circulation rate (F) increases or the 
difference between concentrations £3 and ^ decreases, the COP of the system decreases. 
In the present study, the increase or decrease in F does not provide any information about 
the COP of the system. This is because T g is not constant. For the water-sulfuric acid 
pair the specific circulation rate increases rapidly from 4.341 kg/kg of vapor at T e = 373K 

O 

to a maximum of 110.6 kg/kg of vapor at T g = 338K. The increase in F is due to the 
difference between £3 and as a function of T g . The change in specific circulation rate 
as a function of temperature for both systems is plotted in Figure 8 . 

The maximum cooling mode COP is obtained when Tg is around 353K. The 
cooling mode COP is the ratio of the heat absorbed at the evaporator to that supplied at 
the generator. With an assumed constant load, the heat absorbed at the evaporator (Q e ) 
does not change and, therefore the only variable is the heat supplied to the generator (Q e ) 

O 

Q g depends on several factors. As seen from Equation 8 , h 2 , h 3 , m 3 , and m 2 are the key 

parameters that affect Q g . Q g decreases as I 13 decreases and h 2 increases. As T e 

© 

decreases from 373K to 338K, both 113 and h 2 decrease. The ratio of h 2 to 113 is largest at 
a generator temperature of about 353K; hence the COP is maximum. 

A plot of the heating-mode COP as a function of the generator temperature is 
shown in Figure 9. For the case of water-sulfuric acid, the COP increases from 1.788 at 
Tg = 373K to 1.805 at Tg = 353K, and then decreases. Similarly, the COP of the water- 
lithium bromide system increases steadily from 1.801 to 1.817 for the same temperature 
range. The increase in COP in the former case can be explained by the continuous 
increase in heat of dilution evolved at the absorber, Q a . This heat is largely affected by the 
enthalpy of the weak solution entering the absorber. On the other hand, as the 
temperature lift decreases, Q a increases, and hence, overall COP increases. 
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AHP Performance u sing the Water-Magnesium Chlnririp Pair 

The performance of water-magnesium chloride as an alternative refrigerant- 
absorbent pair is summarized in Tables 11 and 12 . The temperatures for the performance 
calculation are T a =298K, T e =283K, and T C =298K. The generator temperature varies 
from 321K to 315K. The concentration of the strong solution leaving the generator (£ 3 ) is 
greater than that leaving the absorber (^). This limits the range over which T g can be 

varied. The maximum temperature of the generator is limited to 32 IK, since the generator 
and condenser pressures must be equal. 

Table 11 shows the results obtained using water-magnesium chloride as the 
working pair. The performance of this pair is compared with those obtained using water- 
lithium bromide for the same heat pump conditions. The results of the water-lithium 
bromide pair are shown in Table 12. The cooling COP of both pairs are compared in 
Figure 10 . In the cooling mode, the COP of the water-magnesium chloride pair remains 
steady at 0.76 over the entire range. The COP using water-lithium bromide decreases 
slightly from 0.904 to 0.794 for the same range. 

Figure 1 1 shows the comparison of the specific circulation rates. As can be seen, 
the specific circulation rates increase with decreasing temperature lifts. The sudden 
increase in F for water-magnesium chloride is due to the appreciable change in mass flow 
rate at the absorber outlet (mi). The only disadvantage with higher specific circulation 
rates is that more solution must be circulated per unit mass of the refrigerant. 

The heating mode COP of both pairs are shown in Figure 12. The COP for water- 
magnesium chloride remains nearly steady at 1.76, while that for water-lithium bromide 
decreases from 1.904 to 1.794. The water-magnesium chloride COP remains constant 
because Q c , Q a and Q g do not change appreciably. 
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Conclusion 


The AHP analysis has resulted in the discovery of some interesting features of 
water-sulfuric acid and water-magnesium chloride as refrigerant-absorbent pairs, as 
compared to the water-lithium bromide system. These are summarized as follows: 

1. The water-sulfuric acid pair can be used as a satisfactory replacement for water- 
lithium bromide pair in the targeted temperature range. 

2. Unlike the water-magnesium chloride and water-lithium bromide pairs, the water- 
sulfuric acid pair does not undergo crystallization. 

3. The corrosiveness and toxicity of sulfuric acid must be taken into consideration in 
the design of a practical heat pump. 

4. At lower temperatures, the performance of the AHP using water-magnesium 
chloride does not compare well with those obtained using water-lithium bromide. 

5. The operating temperature and concentration ranges for the water-magnesium 
chloride pair are limited by its solubility. 

6. The computer codes developed can be used to analyze the performance of an AHP 
with any refrigerant-absorbent working fluid pair with only slight modifications. 

The results of this study should help future researchers to design more efficient 
absorption heat pumps. 
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Vapor pressure of water-magnesium chloride solution 
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Enthalpy of Water-Magnesium Chloride Solution 
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Enthalpy of Water-Magnesium Chloride Solution 
T(K)=[373-473], m(moI/kg)=[0.1-3.6] 
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= solution temperature (K) [373 < T < 473] 

= molality of electrolyte (mol / kg) [0.1 <m <3. 6] 




Hot Plate 



33 


Absorption Heat Pump Configuration 
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Cooling Mode COP Comparison [H20-MgCl2 Vs. H20-LiBr] 
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Variation in Specific Circulation Rates with Temperature 
[H 2 0-MgCl 2 Vs. H 2 Q-LiBr] 
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Heating Mode COP Comparison [H 2 0-MgCl 2 Vs. H 2 0-LiBr] 
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Design of an ElectroHydrodynamic Pump for the Absorption 

Heat Pump 



Design of an ElectroHvdrodvnamic Pump for the Absorption 


Heat Pump 
Introduction 

Several thermal control techniques have been developed for thermal transport and 
heat rejection in space applications (Kerrebrock, 1986). One method is the absorption- 
cycle heat pump system. It is generally understood that the absorption heat pump offers 
significant benefits in terms of mass and complexity reduction of the radiator area, tighter 
thermal control, and relatively high efficiency. One essential component of the absorption 
heat pump is a highly reliable pump which is needed to circulate the working fluids. 
However, major problems limiting the reliability of currently used pumps include wear in 
the bearings, seals, and in moving parts of the pump. The electrohydrodynamic (EHD) 
pump is an innovative technology which, through proper development, can provide a 
solution to this problem. 

As early as 1937, J. Hartman and F. Lazarus investigated the EHD flow of 
mercury between two flat plates assumed infinite in two dimensions (Samaras, 1962). 
This fluid flow arrangement, known as Hartman-Poiseuille flow, was used in the 
theoretical development of C. Beal (1967) and J. Melcher (1981). Since the work of 
Hartman and Lazarus, several types of EHD pumps have been studied (Baker, 1987), but 
the most appropriate type of EHD pump for aqueous fluids to date is a conduction, or 
direct-current Faraday type. Also called simply an electromagnetic, or "EM" pump, this 
type works on the principle of current flow across the fluid perpendicular to an external 
magnetic field. A resultant body force is exerted on the fluid in a direction mutually 
perpendicular to the current and to the magnetic field in the sense of a left-handed 
coordinate system. 
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Liquid metal EM pumps based on this principle have been available since the early 
1950's and are used primarily in nuclear reactor cooling facilities. Pumps of this type have 
been built for volumetric capacities as high as 3500 mVhour and for heads approaching 
2.4 X 10 5 Pa (Andrev, 1981). 

Investigations are currently being conducted by others on the use of 
electrohydrodynamic pump for sea water thrusters (Doss, 1992). Also, this unique pump 
is under investigation by the U.S. Navy for the development of a silent submarine 
propulsion system (Swallom, Sadovnik, Gibbs, Gurol, Nguyen, van den Bergh & Hugo, 
1991). 

The EHD pump is considered to be a reliable means of circulating the working 
fluid in a heat pump system because it has no moving parts, seals, or bearings. In a pump 
of this type, the pumping force for moving the fluid does not originate from a piston, 
diaphragm, or impeller. Instead, the force is exerted directly on the fluid. The EHD 
pump, therefore, offers great promise for application in an absorption heat pump system. 
Furthermore, the added feature of being able to repair or service these devices without 
breaking the fluid circuit makes this pump more attractive than conventional pumps. 
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Theoretical Modeling of an EHP Pumn 


A continuum model is used to describe the characteristics of the EHD pump, 
according to the following assumptions: 

1. Incompressible, viscous, steady, and fully developed flow. 

2. Uniform electromagnetic field applied in a direction perpendicular to the flow. 
The physical configuration is illustrated in Figure 1. The governing equations for 

the fluid, in normalized form (Melcher, 1981) are as follows: 


V • H = 0 


V x E = - — 


at 


0) 

( 2 ) 


V x H = — CT(E x V x H) 
X 


t j j 

— + V • Vv + VP = -a (E + V x H) x H H — — V 2 V 

X MI ^MI T v 


(3) 

(4) 

'MI “MI <- v 

V • V = 0, (5) 

where, H, E, V, P, and t are electric field intensity, magnetic field intensity, fluid velocity, 
fluid pressure, and time, respectively. The characteristic times: 


T 


1 rl 2 

-> % = 

U rj 

CT o4 0 l 2 , X MI = Vp i2 /p h2 


( 6 ) 

( 7 ) 


are, respectively, the transport time, viscous diffusion time, magnetic diffusion time, and 
the magneto-inertial time. 1 is the characteristic length, U is the typical material velocity, 
p is the fluid density, and rj is the fluid viscosity. It is assumed that the fluid is an ohmic 
conductor with characteristic conductivity of a 0 and essentially the permeability of free 
space, p 0 . 
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In typical laboratory experimentation involving the EHD flow of electrolytes, 
liquid metals, or plasmas, the magnetic diffusion times are short compared to the time of 
interest (e.g., x m «x). Hence, the magnetic Reynolds number (Re m = x m /x) plays 

almost no role in most parts of the flow regime. It is, therefore, appropriate to assume 
that the induced currents on the right-hand side of Equation 3 are negligible. The 
magnetic field is imposed by means of currents in the external windings. With the above 
assumptions, the governing equations are reduced to: 


V x E = 


dH 


dt 

V • J = 0; J = o(E + V x p o H) 
dv 

p(— + (v • V)v) + vp = j x p o H + -nv 2 v 


( 8 ) 

( 9 ) 

( 10 ) 


Locally, the axial variations of electric field and current density are assumed to be 
negligible in comparison with their variations perpendicular to the flow. This assumption 
may not be accurate when strong variations exist in the magnetic and/or velocity flow 
field, or when abrupt changes in the boundary conditions occur, as may happen near the 
pump entrance and exit. However, this analysis is limited to the steady and fully 
developed section of the flow. 

For the task of finding a two-dimensional, steady-flow solution, Equation 10 reduces 
to 


+ _ g( p 0 HJ 2 V y = p o PH 0 E + J. dP 

3x 2 dz 2 r) r\ 2 t| dy ‘ 

The boundary conditions of this flow require that the velocity of the fluid be 


( 11 ) 

zero 


at the fluid/wall interface. The volumetric flow rate Q is calculated from: 


Q = J V,dA c 


( 12 ) 
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where Aq is the cross-sectional area of the pump. 

The electrical power applied to the pump (P e i ec ) is obtained from 

P 8lec = J (J • E)dV, (13 ) 

V 

where V is the active volume of the pump. Since some part of the electrical power is 
consumed in Joule heating, the mechanical power delivered to the fluid is, thus lower than 
the electrical power applied to the pump. The power imparted on the fluid (P fluicl ) can be 
determined from 

P nuid * J(F ' V y )dV, ( 14 ) 

V 

where F is the Lorentz force defined as F = (J x B). The actual power produced by the 
Pump (P out ) is calculated from: 

P out - Q • AP, (15) 

where AP is the differential pressure across the pump. The actual power produced by the 
pump is lower than the power imparted on the fluid due to the fact that some of this 
power is dissipated as internal viscous loss. 

Three different types of efficiency can be defined to provide an insight into the 
behavior of the pump. They are: 

1. Electrical efficiency, defined as. 


^lelec 


f (F • V)dV 

P J 

•Ff luld _ V 

P elec J (J ' E)dV 


V 

2. Internal efficiency, defined as: 


(16) 




p out = Q • AP 
p nuid J ( F • V y ) dV 


(17) 


48 



3. Overall efficiency, defined as: ' 


*lo 


Pout = Q • AP 

P.lec j (J • E)dV 


= "Hint 


"Helec 


These three efficiencies are discussed in detail in later sections. 



Numerical Method 


Referring to Equation 1 1, 
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Vj j is the velocity calculated from Gauss-Seidel iteration at node i, j 
v i,j is th e approximation of V y by the SOR method at node i, j. 

<o is a relaxation parameter. 

<k> is the component of the kth iteration. 

The optimal value of the relaxation parameter, co, is found empirically by 
performing calculations with different values of co for coarse grids until the optimal value 
of a is obtained. With the optimal co, the corresponding optimal value of a for finer 
grids can then be found by applying the following formula (Strikwerda, 1989) : 

2 


where C is a constant obtained experimentally, and h is a distance between iterative nodes. 

A computer program (Appendix B) was written to implement the iterative scheme. 
Numerical results from this program are then transferred to a graphics software package 
for graphical presentation. 
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Experimental Apparatus and Instrumentation 


Experiments were performed to verify the theoretical results obtained from 
modeling the electromagnetic flow. A copper sulfate solution was used as the working 
fluid for model verification, with properties as shown in Appendix B. 

Experimental Apparatus 

The experimental apparatus, illustrated in Figure 3, consists of a reservoir, 
rectangular conduit, EHD pump, manometer, rotameter, and secondary centrifugal pump. 

As can be seen from Figure 3, the rectangular conduit carries the fluid from the 
reservoir to the EHD pump, and is connected by the return conduit (7.875 mm ID). The 
flow rate was monitored by a rotameter connected to the end of the return conduit. A 
valve at the rotameter was used to control the flow rate of the system. The secondary 
centrifugal pump was connected from the reservoir to the rectangular conduit in parallel to 
the direct connection between the reservoir and the rectangular duct. This secondary 
pump was used to drain air bubbles from the system during the startup process. Startup 
and operational procedures are discussed later in this chapter. 

The rectangular conduit, with internal dimensions of 6.35 mm by 25.4 mm, was 
made of Plexiglas. Plexiglas was chosen because, in addition to being a good electrical 
insulator, it also allows visual observation of the pumping process. The various Plexiglas 
parts were assembled using dichloroethane solvent. 

Two 6.35 mm by 101.60 mm slots were provided on each side of the conduit at 
the pump section to accommodate the copper electrodes, which were affixed to the 
Plexiglas conduit with epoxy. 

Figure 4 illustrates the EHD pump used for the experiments. The EHD pump 
envelope is by definition the volume between the electrodes and the magnetic poles. Some 


52 



consideration was given to the fact that the theoretical modeling was based on the 
electromagnetic flow in the folly developed region. The pump section, thus, should be as 
long as possible to comply with the theoretical model. The length of the pump section 
was chosen to be 101.6 mm (four times greater than the width of the channel). 

The copper electrodes, shown in Figure 4, were connected to the power supply via 
three #18 AWG insulated wires to ensure even distribution of the electric field across the 
pump section. The surfaces of the electrodes were sanded and polished before installation 
in the conduit. A 15 ampere, 0-50 volt DC power supply (Accurate Electronics Co., 
model PP-1725) was used with the electrodes. 

The magnetic circuit was designed to provide a maximum of 0.4 teslar magnetic 
flux density in the vertical direction across the pump section. The geometry of the steel 
core and the arrangement of the coil winding were configured to minimize magnetic 

flinging outside the pump. Figure 5 illustrates the details of the electromagnet used in the 
test setup. 

Since a larger wire requires higher current at lower voltage and vice versa, an 
analysis (Appendix C) was performed to select the proper wire size for operation at the 
voltage and current rating of the available power supply. According to the analysis, #18 
AWG wire with 840 turns is sufficient for this application. 

A 10 ampere, 0-55 volt DC power supply (Epsco Incorporated., model PS-5) was 
used for the electromagnet power supply. The actual magnetic flux density was measured 
using a gaussmeter (General Electric, model 416X37). This reading was then used as a 
standard for the entire set of experiments. Comparison of the calculated and measured 
magnetic flux density is given in Appendix D. 

A manometry system was used to measured the difference between the suction and 
discharge heads. Two static pressure taps of 2.23 mm diameter at the entrance and exit of 
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the EHD pump were connected to the manometer tubes (glass tubes of 4.5 mm internal 
diameter) by 3. 175 mm internal diameter flexible plastic tubing. 

The flow measurements at high flow rate were acquired using a Dyna VFB-141- 
BV, 0.5-5 gpm (31.5-10- 6 -315.410- 6 m 3 /s) rotameter. The rotameter was calibrated for 
the test fluid (Appendix E). To obtain accurate readings at flow rates lower than 3 1 .5- 10’ 

6 m 3 /s, the flow measurements were made by measuring the total volume collected in a 
given time interval. 

Initially, a pitot tube was installed at the end of the pump to measure the local 
velocity of the fluid. However, this method was abandoned due to unsatisfactory 
differential pressure measurement. 

Two multimeters (Micronta Electronics Inc., model 22- 185 A) were used to 
measure the voltage and current applied to the electromagnet and electrodes before each 
experiment was performed. Once the experiments were started, the voltage and current 

readings were then monitored using the voltmeters and ammeters provided with the power 
supplies. 


Experimental Procedure 

After apparatus assembly, leak-checking, and instrumentation verification, the 
experimental system was certified "test-ready," and all was in order to begin the test runs. 
Before beginning, however, some precautions were taken to assure the integrity of the test 
data. These are summarized as follows: 

1. The valve at the rotameter must be fully opened when the centrifugal pump is 

turned on in order to prevent hydraulic pressure from building up and rupturing the 
conduit. 

2. Prolonged use of the electromagnet may cause a damage to the magnetic circuit 
due to high temperatures generated in the interior coils. Thus, coil temperature was 


54 



gauged by noting the rising resistance of the magnetic circuit. In this way the resistance of 
the magnetic circuit was monitored continuously and maintained lower than 6 ohm during 
the experiments. 

3. The resistance across the electrodes must be regularly monitored. If the 
resistance is higher than the initial value (8.8 ohm), the electrodes must be cleaned. This 
increase in resistance cross the electrodes is due to the uneven deposition of copper and 
impurities on the cathode. 

After having incorporated the above precautions into the test operation, the 
following experimental procedure was implemented: 

1. The test loop was filled with the test solution. This was done by pouring the 
test solution into the reservoir until the level of the test solution was approximately 10 cm 
above the conduit level. The valve at the rotameter was then slowly opened to allow the 
solution to flow into the conduit. Trapped air inside the conduit was drained by closing 
the ball valve located near the reservoir and turning on the centrifugal pump to help 
circulation and to force air to exit the system at the reservoir. 

2. The magnetic circuit was energized. Multimeters were used to monitor the 
current and voltage applied to the electromagnet. The supply source was adjusted until 
the power necessary to provide the desired flux density was obtained. The power supply 
was then turned off. 

3. The power supply to the electrodes was turned on. The power supply was 
adjusted until the desired voltage across the electrodes was obtained. Two multimeters 
were used to monitor the voltage and current supplied to the electrodes during this 
process. The power supplied was then turned off after the desired power was set. 

4. For each position of the valve settings at the manometer (from fully-closed to 
fully-opened), the magnetic field and the electric field were simultaneously energized. The 
ball valve at the reservoir was then opened. The differential pressure head across the EHD 
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pump was measured by the manometer. The lower surfaces of the suction pressure and 
the discharge pressure meniscuses formed in the manometer tubes were marked. The 
height difference was then recorded. The flow rate was also observed directly at the 
rotameter for high flow rates. For lower flow rates, where the rotameter reading was 
inaccurate, the time interval required to fill a calibrated container was recorded, and the 
flow rate was calculated from this data. 

5. The experiment was performed for various magnetic flux densities and voltages 
applied across the electrodes. Table 1 summarizes the various conditions of magnetic flux 
density and voltage across the electrodes used in the test runs. 
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Results .and Discussion 


Numerical results 

The following sections discuss the numerical results obtained for the velocity 
profile, pressure-flow rate relationship, and efficiency of the EHD pump as a function of 
magnetic flux density and conduction current. 

1. Numerical Results for the Velocity Profile 

The numerical methods developed earlier have been used to determine the flow 
characteristics of the EHD pump. The numbers of nodes in the x and z directions were 
initially chosen to be 11 and 41, respectively. The dependence of the accuracy of the 
results on the number of nodes chosen was then investigated. It was found that the results 
varied by less than 1% when a 21 X 81 -node configuration was used, and more than 3% 
for a 5 X 17-node configuration. Hence, the 1 1 X 41 -node configuration was retained for 
this analysis. 

The relaxation factor ( a ) is an empirically determined quantity of the calculational 
process. Since the solution converges rapidly, the relaxation factor can be found directly 
without using the scaled-down configuration (lesser number of nodes). It was found that 
the most suitable value of the relaxation factor for the nodal configuration used (11 X 41) 
was 1.7. 

Since the mathematical model has no closed-form solution, verification of the 
results was accomplished by modifying the two-dimensional finite-difference scheme 
developed earlier to a one-dimensional model, and then comparing the results obtained 
with the analytical closed-form solution. Comparison of the numerical and closed-form 
solutions indicated an error of less than 1%. 
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An example of the numerical results obtained for the axial, fully-developed velocity 
profile is shown in Figure 6. The vertical plane indicates the axial component of the 
velocity (Vy). 

The Lorentz force in the axial direction (J z X B) acts as the driving force for the 
pump. This force, however, is not uniform over the entire cross-section of the pump flow 
area, due to the non-uniform distribution of the current density. It can be seen from 
Equation 7 that the local current density depends upon the local fluid velocity. For an 
EHD pump, a higher fluid velocity corresponds to a lower current density. The force 
acting on the fluid near the walls, where the local velocity is lower, is higher than that at 
the core of the duct. This increase in the force causes a rapid change in the velocity 
profile near the walls. 

From the above discussion, it appears that two factors, namely, the voltage across 
the electrodes and the magnetic flux density, are the parameters controlling the 
characterization of EHD pump performance. Figure 7 (a,b and c) illustrates the resultant 
velocity profiles for various applied voltages with constant magnetic flux density and 
pressure. The voltage across the electrodes exhibits a synergistic effect on the velocity 
profiles as indicated in the figures for voltage variations from 10 V to 20 V. In Figure 8 
(a, b and c) , the effect of varying the magnetic flux density is shown. As with the 
previous case, the velocity profiles intensify as the flux density increases from 0.2 T to 0.4 
T. 

2. Numerical Results for Pressure-Flow Rate Performance Characteristics 

The pressure-flow rate relationship is also of primary interest in characterizing the 
performance of the pump. A parametric study was performed by varying the differential 
pressure while the pump power input was held constant. The analysis indicates that the 
theoretical pressure-flow rate relationship for an EHD pump is similar to that of a 
conventional centrifugal pump. 
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Figures 9 illustrates the pressure-flow rate performance characteristic for various 
magnetic flux densities with constant voltage maintained across the electrodes. As 
indicated in the figure, an increase in the flux density increases the pump output. 

The effect of the voltage applied across the electrodes on the pressure-flow rate 
performance characteristic is illustrated in Figure 10. As can be seen from this figure, the 
pump output increases with an increase in the applied voltage. However, it should be 
noted that the physical properties of the fluid limit the voltage across the electrodes. This 
limitation is discussed later in this chapter. 

3. Numerical Results for Pump Efficiency 

Since efficiency is a parameter of primary importance in characterizing the 
performance of an EHD pump, the efficiency was determined using two different 
methods. 

In the first approach, the effect of magnetic flux density and voltage applied to the 
electrodes on pump efficiency was investigated while maintaining a constant flow rate. 
This method is useful in determining the optimal input condition for a given flow rate. 

In the second method, the effects of flow rate on pump efficiency were studied for 
a given voltage and magnetic flux input. This method provides a convenient way of 
comparing the numerical results with the experimental results. The comparison of the 
numerical and experimental results are discussed later in this chapter. 

<• Effects of Magnetic Flux Density and Voltage Across the Electrodes on EHD 
Pump Efficiency . 

Three different types of efficiencies have been previously discussed. Figure 1 1 
illustrates the effects of the magnetic flux density and the voltage across the electrodes on 
the internal efficiency of the EHD pump for a constant volumetric flow rate. As can be 
seen from Figure 11, an increase in either the magnetic flux density or the voltage across 
the electrodes increases the internal efficiency of the pump. This is because, at constant 
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flow rate, the viscous loss in the fluid is almost constant. For a higher magnetic flux 
density or an increase in the voltage across the electrodes, the pump power input increases 
while the losses remain nearly constant. This improves the internal efficiency of the EHD 
pump. 

The effects of varying magnetic flux density and voltage across the electrodes on 
the electrical efficiency are shown in Figure 12. As indicated in the Figure, high electrical 
efficiency can be obtained using higher magnetic flux density and lower voltage across the 
electrodes. This phenomenon can be explained by considering local electrical efficiency 
given as: 


Klee ] 


B • V. 


B • V„ 


local 


( 17 ) 


E v/w 

where B is the magnetic flux density; Vy the local velocity of the fluid and E the electric 
field strength. E is again defined as the ratio of the voltage across the electrodes to the 
width of the pump (E = v/w). For a given flow rate, the local fluid velocity is almost 
constant. The local electrical efficiency thus depends on the ratio of the flux density to the 
voltage across the electrodes. 

It should be noted that the local electrical efficiency is not uniform over the entire 
cross-sectional area of the pump. Equation 17 indicates that the local electrical efficiency 
depends on the local fluid velocity. The local electrical efficiency of the fluid at the core 
is, therefore, higher than that near the walls. Figure 13 illustrates the overall efficiency of 
the pump. The overall efficiency of the pump is the product of the integrated electrical 
and internal efficiencies. An increase in the magnetic flux density will increase both the 
electrical and internal efficiencies, therefore increasing the overall efficiency of the EHD 
pump. An increase in the voltage across the electrodes will increase the internal efficiency 
but decrease the electrical efficiency. The maximum overall efficiency of the pump is 
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therefore dependent on the particular pump configuration, namely, flow rate, voltage 
across the electrodes, and magnetic flux density, 
ii. Effect of Flow Rate on EHD Pump Efficiency 

The effect of volumetric flow rate on the efficiency of the pump is another 
important characteristic which must be understood to properly characterize pump 
performance. A parametric study was conducted by varying both the volumetric flow rate 
and the pump power input. Figures 14 (a and b) illustrate the relationship between 
internal efficiency and volumetric flow rate for various input conditions. As can be seen 
from the figure, the internal efficiency decreases linearly from its maximum (no-flow 
condition) to zero (maximum flow rate). This is due to the fact that the viscous loss of the 
fluid is directly proportional to the volumetric flow rate. For the maximum flow rate, all 
of the mechanical power delivered to the fluid is used to overcome the viscous loss. 
Therefore, no useful output is produced. The internal efficiency of the pump is thus zero 
for the maximum delivery condition. Since the viscous loss converges to zero as the flow 
rate approaches the no-flow condition, the internal efficiency is at its maximum. 

Figure 15 (a and b) shows the electrical efficiency-volumetric flow rate relationship 
for the EHD pump. Increasing the flow rate results in an increase in the mechanical power 
delivered to the fluid but does not change the input power to the system. Therefore, the 
electrical efficiency increases as the flow rate increases. 

In Figure 15a, it is important to note that as the magnetic flux density increases, 
the slopes of the electrical efficiency-flow rate curves (each at a constant B) increase. But 
as seen from Figure 1 5b, the slopes of the electrical efficiency-flow rate curves (each at a 
constant v) decrease as the applied voltage across the electrodes increases. 

Overall efficiency curves for the EHD pump are shown in Figures 16 (a and b). 
The overall efficiency as defined earlier is the product of the internal and electrical 
efficiencies. Figures 16a shows that for a given voltage across the electrodes, the overall 
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efficiency increases from zero to a maximum and then decreases to zero as flow rate 
increases from the no-flow condition to the maximum-flow condition. As the flux density 
increases, the maximum value of overall efficiency also increases. It should be noted that 
the maximum values of the overall efficiency occur at higher flow rates when the flux 
density increases. 

Figure 16b shows the theoretical relationship between overall efficiency and flow 
rate for constant magnetic flux density and varying voltage. For a given value of voltage 
applied across the electrodes, the overall efficiency increases form zero to its maximum 
and then decreases to zero as the flow rate increases. As the voltage increases, the 
maximum value of the overall efficiency obtained remains approximately the same but 
occurs at higher flow rate values. 


Experimental Results 

Since repeated attempts to measure the velocity profile yielded unsatisfactory 
results, the following sections will focus only on the experimental results obtained for the 
pressure-flow rate relationship and efficiency determinations for the EHD pump. 

1. Experimental Results of Pressure-Flow Rate Relationship 

The effects of the magnetic flux density and voltage across the electrodes have 
been discussed in previous sections. Figures 17 to 28 show the comparison between the 
experimental and theoretical pressure-flow rate relationship for various input conditions. 
The various input conditions are summarized in Table 2. 

As seen from Figures 17 to 28, the experimental values agree well with the 
predictions for low flow rates. At higher flow rates, the experimentally measured values 
were lower than the predicted values. The aberrations observed are due to the assumption 
that the flow is fully developed. Since the fluid at low flow rate tends to become fully 


62 



developed faster than the fluid at high flow rate, the theoretically predicted values are 
more accurate at the low flow rates than at the high flow rates. 

The difference in the measured and predicted pressure at high flow rates may also 
be due to the non-linear losses in the system. These non-linear losses are the those caused 
by local turbulence. 

For higher values of applied voltage across the electrodes, the efficiencies obtained 
experimentally were lower than the predicted values. The variations observed were due 
to increasing current fringing at higher voltages. This resulted in the decrease of the 
useful conduction current in the system. 

Physical properties of the working fluid also limit the applied voltage. High 
voltage applied across the electrodes caused an uneven deposition of copper on the 
cathode, thus increasing the electrode surface resistance. This effect degraded the pump 
power output rapidly. Since this phenomenon is pronounced for applied voltages higher 
than 20 V, no experiments were performed in this voltage range. 

As the magnetic flux density increased, it was expected that the fringing effect 
would cause variations in the experimental results, in the same manner that current 
fringing did. However, comparison of experimental and theoretical data does not 
conclusively substantiate this theory. 

2. Experimental Results for Overall Efficiency 

Since internal and electrical efficiency could not be directly measured, the overall 
efficiency was the most important parameter of interest. The theoretical effects that the 
magnetic flux density, the voltage across the electrodes, and the volumetric flow rate have 
on the overall efficiency have been discussed previously. 

The experimental and theoretical values of overall efficiency are shown in Figures 
29 to 40 for the various input conditions. The input conditions are listed in Table 3. 
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As can be seen from Figures 29 to 40, the experimentally obtained efficiencies 
agree well with the predicted values at low flow rates. However, the experimental results 
were lower than the numerical results for higher flow rates. This is due to the fact that the 
actual differential pressure across the pump is lower than the theoretical value, when the 
flow rates become higher. 

The current fringing effect could also cause the actual overall efficiency to be 
lower than the theoretically predicted value. As discussed earlier, higher voltage across 
the electrodes increases the current fringing, thus reducing the power output of the pump. 
Although theoretically, the magnetic fringing should decrease the pump performance, 
experimental data does not substantiate this fact. 
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Conclusions and Recommendations 


An investigation of the operational characteristics of an electrohydrodynamic pump 
was performed in this study. The numerical and experimental results obtained provide an 
insight into the electrohydrodynamic characteristics of the fluid inside the pump. The 
results obtained are summarized as follows: 

1. The numerical mathematical model developed based on continuum 
electromechanics provides reasonable results within the limitations of the assumptions 
made. 

2. The operational characteristics of the EHD pump as predicted by the model 
were found to be similar to those of a conventional centrifugal pump. 

3. The predicted overall behavior of the EHD pump is comparable to that of 
turbomachinery hardware; however, the overall efficiency of the electromagnetic pump is, 
as predicted, much lower. 

4. The pressure-flow rate performance characteristics and overall efficiency of the 
pump obtained experimentally agree well with the theoretical results. 

Although the model developed during this study yielded satisfactory results, better 
results could be obtained with a more sophisticated model based on a three-dimensional 
mathematical model. 

The small size of the experimental apparatus used in this study caused some 
inherent limitations in the accuracy of the measured quantities. In the future, these 
limitations could be eliminated by utilizing a pump of larger capacity. 

The electrohydrodynamic pump is very sensitive to fluid properties. Because of 
this, care must be exercised in selecting the feasible range of operational parameters and 
appropriate electrode types for a given application. As a recommendation for future 
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research, it is proposed that additional effort be expended to explore the variety of fluids 
available, with appropriate electrode combinations, for use in an EHD pump. 

Also, as the overall efficiency of the pump is largely dependent on the magnetic 
flux density, superconductors could be used to improve pump performance. 

The development of the EHD pump for circulating the working fluid in the 
absorption heat pump system can be considered a high-risk, high-payoff application of 
advanced technology to provide revolutionary improvement in reliability and 
maintainability of the heat pump system. The results of the present effort should help 
researchers understand the fundamental performance characteristics of this type of pump 
in order to develop an improved electrohydrodynamic pump for future applications. 
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EHD Pump Layout used in the Experiments 




Details of the Electromagnet 
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An Example of the Velocity Profile Obtained Numerically 
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The Velocity Profile for V = 10 V (B = 0.2 T, P = 30 Pa) 
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The Velocity Profile for V = 15 V (B = 0.2 T, P = 30 Pa) 
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The Velocity Profile for V = 20 V (B = 0.2 T, P = 30 Pa) 
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Figure 10 

The Velocity Profile for B = 0.2 T (V = 15 V, P = 30 Pa) 
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The Velocity Profile for B = 0.3 T (V - 15 V, P - 30 Pa) 
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The Velocity Profile for B = 0.4 T (V = 15 V, P = 30 Pa) 




Figure 13 

The Pressure - Flow Rate Relationship at Constant Voltage (V = 15 V) 




001 



82 


The Pressure - Flow Rate Relationship for Constant Magnetic Flux 
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Contour Plot of the Internal Efficiency at Constant Volumetric Flow Rate (Q = 1 x 1(H m 3 /s) 
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Contour Plot of the Overall Efficiency at Constant Vol 
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Figure 19 

Internal Efficiency as A Function of Magnetic Flux Density at Constant Voltage (V — 15 V) 
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Electrical Efficiency as A Function of Magnetic Flux Density at Constant Voltage (V - 15 V) 
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Overall Efficiency as a Function of Magnetic Flux Density at Constant Voltage (V = 15 V) 
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Figure 23 

Overall Efficiency as a Function of Applied Voltage at Constant Magnetic Flux Density (B 
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Pressure - Flow Rate Relationship Obtained Numerically and Experimentally for V = 10 V and B = 0.1 T 
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Figure 25 

Pressure - Row Rate Relationship Obtained Numerically and Experimentally for V = 10 V and B = 0.2 T 
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Figure 26 

Pressure - Flow Rate Relationship Obtained Numerically and Experimentally for V = 10 V and B = 0.3 T 
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Pressure - Flow Rate Relationship Obtained Numerically and Experimentally for V = 10 V and B = 0.4 T 
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Pressure - Flow Rate Relationship Obtained Numerically and Experimentally for V = 15 V and B = 0.1 T 
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Figure 29 

Pressure - Flow Rate Relationship Obtained Numerically and Experimentally for V = 15 V and B = 0.2 T 
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Figure 30 

Pressure - Flow Rate Relationship Obtained Numerically and Experimentally for V = 15 V and B 
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Pressure - Flow Rate Relationship Obtained Numerically and Experimentally for V = 15 V and B - 0.4 T 
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Figure 32 

Pressure - How Rule Rela.ionship Obtained Numerically aud Eiperimeutally for V = 20 V and B - 0.1 T 






Figure 34 

Pressure - Flow Rate R.| a ,i„„ s hip Obtained Numerically and Experimentally for V = 20 V and B = 0.3 T 
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Figure 35 

Pressure - Flow Rate Relationship Obtained Numerically and Experimentally for V = 20 V and B = 0.4 T 
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Overall Efficiency Obtained Numerically and Experimentally for V = 10 V and B = 0.2 T 
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Figure 38 

Overall Efficiency Obtained Numerically and Experimentally for V = 10 V and B = 0.3 T 
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Overall Efficiency Obtained Numerically and Experimentally for V — 10 V and B = 0.4 T 
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Figure 41 

Overall Efficiency Obtained Numerically and Experimentally for V = 15 V and B = 0.2 T 
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Overall Efficiency Obtained Numerically and Experimentally for V = 15 V and B = 0.3 T 
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Overall Efficiency Obtained Numerically and Experimentally for V = 15 V and B = 0.4 T 



Experimental Results 


1 


T 


i 

I 


in 

o 

o 

o 


□ 


L. 

* 

o 

o 

o 




o 


o 

in 


o 


CO 


ro 

£ 


I o 


X 

O 


114 


I 


Figure 44 

Overall Efficiency Obtained Numerically and Experimentally for V = 20 V and B = 0.1 T 






Figure 45 

Overall Efficiency Obtained Numerically and Experimentally for V = 20 V and B = 0.2 T 
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Figure 46 

Overall Efficiency Obtained Numerically and Experimentally for V = 20 V and B = 0.3 T 




Overall EfTiciency Obtained Numerically and Experimentally for V = 20 V and B = 0.4 T 
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Appendix A 


Computer Program 

'FINITE DIFFERENCE FOR EHD PUMP 
'WRITTEN BY WORANUN WORAMONTRI 

CLEAR 
SCREEN 0 
COLOR 4 , 7 
CLS 

LOCATE 3 , 2 0 

PRINT "Finite difference for EHD Conduction Pump" 
LOCATE 5, 37 
PRINT "By" 

LOCATE 7 , 30 

PRINT "Woranun Woramontri" 

'MENU 

COLOR 1, 7 
LOCATE 11, 20 

PRINT "1. Find Velocity Profile Only" 

LOCATE 13, 20 

PRINT "2. Find P-Q Relationship Only" 

LOCATE 15, 20 

PRINT "3. Find Velocity Profile and P-Q Relationship" 
10 LOCATE 20, 20 
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COLOR 6 , 7 

INPUT "Please select (1-3) ", menu 
COLOR 1, 7 


SELECT CASE menu 


CASE 1 

CLS 

LOCATE 5, 24 
COLOR 3, 8 

PRINT " VELOCITY PROFILE CALCULATION " 

CASE 2 

CLS 

LOCATE 5, 24 
COLOR 12 , 1 

PRINT " PRESSURE- FLOW RATE CALCULATION " 

CASE 3 

CLS 

COLOR 14 , 3 
LOCATE 5, 29 

PRINT " VELOCITY PROFILE " 

LOCATE 7, 37 
PRINT " AND " 

LOCATE 9, 24 

PRINT " PRESSURE -FLOW RATE CALCULATION " 
CASE ELSE 
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GOTO 10 


END SELECT 
COLOR 1, 7 

' INPUT DATA 
LOCATE 12, 20 

INPUT "Number of Vertical Nodes"; nv 
IF nv = 0 THEN nv = 11 
LOCATE 13, 20 

INPUT "Number of Horizontal Nodes"; nh 
IF nh = 0 THEN nh = 41 
LOCATE 14, 20 

INPUT "Distance Between Nodes"; h 
IF h = 0 THEN h = 0.00006392 
LOCATE 15, 20 

INPUT "Length of the Pump"; length 
IF lenth = 0 THEN length = 0.1016 
LOCATE 16, 20 

INPUT "Magnetic Flux Density" ; b 
LOCATE 17, 20 

INPUT "Voltage apply Across the Electrodes"; volt 
SELECT CASE menu 
CASE 1 
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LOCATE 18, 20 

INPUT "Pressure Rise Across the Pump"; dp 
LOCATE 19, 20 

INPUT "File to Write Velocity Profile"; vel$ 
vel$ = vel$ + " .DTA" 

OPEN vel$ FOR OUTPUT AS #1 
dta = l 


CASE 2 

LOCATE 18, 20 

INPUT "Maximum Pressure Rise Across the 
Pump " ; max dp 
LOCATE 19, 20 

INPUT "Data Points Required (for P-Q) " ; dta 
LOCATE 20, 20 

INPUT "File to Write Pressure-Flow Rate"; 

pq$ 

pq$ = pq$ + ".DTA" 

OPEN pq$ FOR OUTPUT AS #2 


CASE 3 

LOCATE 18, 20 

INPUT "Maximum Pressure Rise Across the 
Pump " ; maxdp 
LOCATE 19, 20 

INPUT "Data Points Required (for P-Q)"; dta 
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LOCATE 20, 20 

INPUT "File to Write Pressure-Flow Rate",-pq$ 
OPEN pq$ FOR OUTPUT AS #2 
LOCATE 21, 20 

INPUT "File to Write Velocity Profile"; vel$ 
CASE ELSE 

CLS : PRINT "There is no such a mode!" 

END SELECT 


'INITIALIZE CONSTANTS 

'S.O.R CONSTANTS 
w = 1.7 
ww = 1 - w 


'OPTIMAL S.O.R VALUE 


'ELECTROMAGNETIC/FLUIDIC CONSTANTS 
sig = 4.44 


mu = .000935 

E = volt / ( (nh - 1) * h) 

'FINITE DIFFERENCE CONSTANTS 
disv = (nv - l) / 10 
dish = (nh - l) / 10 

'SYSTEM CONSTANTS 


'FLUID CONDUCTIVITY 
'FLUID DYNAMIC VISCOSITY 
'ELECTRIC FIELD 

'FOR DISPLAY PURPOSE 
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itr = l 


'BEGIN WITH ITERATION 


’VELOCITY PROFILE ARRAY 
DIM u(nh, nv) 

'BOUNDARY CONDITIONS 
FOR i = 0 TO nh - l 
u(i, 0) = 0 
u ( i , nv - l ) = o 
NEXT i 

FOR j = 0 TO nv - l 
u(0 f j) = 0 
u(nh - l, j) = o 

NEXT j 


*********** PHASE I STARTS HERE ********** 
'BEGIN ITERATION 
30 'CONTINUE 

'PRESSURE CONSTANT 
IF menu = 1 THEN 

P = dp / length 

ELSE 

pincre = maxdp / (length 

1 ) ) 

p = pincre * k 


(dta - 
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END IF 


kl 

= (b 

* 2 * 

sig) / mu 

'kl 

k2 

= (P 

/ mu) 

- ( (sig * E * b) / mu) 

' k2 

cl 

= 4 * 

■ kl * 

h * h 

'cl 

c2 

= k2 

* h * 

h 

' c2 


' FILE NAME FOR VELOCITY PROFILE CASE 3 
IF menu = 3 THEN 

decl = INT (k / 10) 

dec2 = 10 * ( (k / 10) - INT (k / 10)) 

FOR i = 0 TO 9 

IF decl = i THEN LET decl$ = CHR$(48 

+ i) 

IF dec2 = i THEN LET dec2$ = CHR$(48 

+ i) 

NEXT i 

ve$ = vel$ + decl$ + dec2$ + " .DTA" 

OPEN ve$ FOR OUTPUT AS #1 

END IF 
CLS 

LOCATE 2, 10 
COLOR 4 , 7 

PRINT "Phase I on progress" 

COLOR 0, 7 
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LOCATE 3, 10 

PRINT "Pressure Rise = "; p * length 
IF menu <> 1 THEN 

LOCATE 3, 30 

PRINT "( Max maxdp; Step = pincre 

* length;")" 

END IF 
LOCATE 4, 10 

PRINT "Calculating Velocity Profile" 

COLOR 6 , 7 
LOCATE 22, 10 

PRINT "Press <s> to stop the iteration" 

COLOR 1, 7 
DIFF = 1 

WHILE INKEY$ <> "S" AND DIFF >= .0000001 AND itr < 
300 

FOR j = 1 TO nv - 2 
FOR i = 1 TO nh - 2 
' Check Convergence 

IF I NT (itr / 10) = itr / 10 AND i = (nh - 1) / 2 AND 
j = (nv - 1) / 2 THEN 

DIFF = ABS (uref - u(i, j)) 
uref = u(i, j) 

END IF 

UU (i , j ) = u (i , j ) 
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•PHASE 


u (i, 

j) = 1 / cl . * (u(i 

+ 1, 

j) + u (i - 1 


u (i, j + l) + 

u (i. 

j - 1) - c2 ) 

u(i, 

j) = ww * UU + w * 

u (i, 

j) 

NEXT i 




NEXT j 





LOCATE 6, 10 
COLOR 6, 7 

PRINT "This is the itr; "th Iteration" 

COLOR 1 , 7 
itr = itr + l 
FOR j = 0 TO 10 
LOCATE 8 + j , 10 

PRINT "value of u ( " ; (nh - 1) / 2; j * disv; ") 

is equal to" 

COLOR 4, 7 
LOCATE 8 + j , 45 

PRINT USING u ( (nh - 1) / 2, j * 

disv) 

COLOR 1 , 7 
NEXT j 
WEND 

I COMPLETE 
BEEP 

LOCATE 2, 10 

PRINT "Phase I completed " 
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IF menu = 2 THEN GOTO 20 ’SKIP THIS PROCESS 

LOCATE 4, 10 

PRINT " „ 

LOCATE 22, 10 

PRINT "Writing file to disk » 

LOCATE 23, 10 

PRINT "Please wait..." 

'WRITING FILE FOR VELOCITY PROFILE 
FOR j = 0 TO nv - 1 
FOR i = 0 TO nh - 1 

PRINT # 1 , USING "##.######*r AAA "; u (i , j) 

NEXT i 
NEXT j 
CLOSE #1 

*********** PHASE II STARTS HERE ********** 


20 'CONTINUE 
CLS 

LOCATE 2, 10 

PRINT "Phase II on progress" 

LOCATE 4, 10 

PRINT "Calculating Volume Flow Rate" 
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sigma = 0 

FOR j = 1 TO nv - 2 
FOR i = 1 TO nh - 2 

sigma = sigma + u(i, j) 

NEXT i 
NEXT j 

q = sigma * h * h 

area = h * h * (nh - 1) * (nv - 1) 

vavg = q / area 

CLS 

COLOR 4, 7 
LOCATE 2, 10 

PRINT "Phase II completed" 

LOCATE 9, 12 
COLOR 1 , 1 

PRINT "Volume flow rate equal to " ; q; " M A 3/S" 

LOCATE 11, 12 

PRINT "Cross-sectional Area equal to area,- "M^2" 
LOCATE 13, 12 

PRINT "Avg. Vel . equal to " ; vavg; "M/S" 

IF menu <> 1 THEN 

PRINT #2, USING »## . ####### AAAA " ; p * 
length; q 

END IF 
itr = 1 
k = k + 1 


129 



IF k < dta THEN GOTO 30 
IF menu = 1 THEN END 
CLOSE #2 
END 


l 


I 

I 
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Appendix B 

Fluid Properties i 

A copper sulfate solution was used as the working fluid for the model verification. 
Its concentration was 133.4 g (CUSO 4 . 5 H 2 O) per liter of water at STP conditions. Its 
other properties, also at STP, are shown in Table B.l. 


Table B. 1 

Properties of Copper Sulfate Solution 


kg 

Density ( — f ) 

m 

1073 

kg 

Dynamic Viscosity ( =— ) 

m * s 

9.35X 10' 4 

„ ^ . . C 2 • S 

Conductivity ( ) 

kg • m 

4 .44 
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Appendix C 

Wire Size Calculation for the Electromagnet 


An analysis was performed to determine the appropriate wire size required for the 
electromagnet. The magnetomotive force required to produce the magnetic flux density 
of 0.4 T was initially calculated from which the appropriate wire size was determined. 
Figure C.l illustrates the geometry of the electromagnet. Table C.l summarizes the 
electrical properties of the solid copper wire considered for this application. 


The magnetomotive force, expressed in terms of ampere-turns, is determined from: 


ni-ba^r*. 


(C.l) 


where N is the total number of turns, I is the current, B is the magnetic flux density, ^ 
is the core cross-sectional area and R (ot is the total reluctance of the magnetic circuit. 

From the geometry of the electromagnet shown in Figure C.l, can be 
expressed as, 


R„ = -^- + 2R, + R, 


(C.2) 


where: 


A co re /f 


R 2 = 


core / s 

l 


Acore/'s 

1, 


R 3 = “ 


A^e^ 
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(C.3) 


The values of lengths 1 1; 1 2 and 1 3 are shown in Figure C.l. The permeability of the 
steel core p s is approximated as: 

H, = 250 

where p 0 is the permeability of free space. 

According to this analysis, 4206 ampere-turns of the magnetomotive force 
required to produce 0.4 T of magnetic flux density at the air gap. 

In order to obtain the required current for each wire size, the number of winding 
turns for each wire size (N) was calculated. The approximation of the winding turns 
be determined from: 


is 


can 


N = 



(C 4) 


where A^, is the total cross-sectional area of the coil and A wirc is the cross-sectional area 
of the wire. The currents required for each wire size were then calculated from: 


I _ NI _ 4206 

N " N (C.5) 

An approximation for the length of the wire used for each case (L coil ) was 
calculated by: 


^ooil N 


(C.6) 


where C^ is the mean circumference of the coil. The total resistance of the coil (R^,) 

was then calculated from the total length and the wire property given in Table C.l. From 

the current and resistance of each wire size, the required voltage is calculated from: 

V =1 r 

V ooil A coil 1 (C 


Figure C.2 shows the required voltages for various wire sizes. The required 
currents for various wire sizes are shown in Figure C.3. From Figures C.2 and C.3, #18 
AWG was considered to be the appropriate wire size for this particular application. 
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Geometry of the Electromagnet 




Properties of Standard Annealed Copper Wire 


l 
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Required Voltage as a Function of Wire Size 
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Required Current as a Function of Wire Size 


Appendix D 

Calculate d and Measured Magnetic Flux Density 


T 


T 



I (Ampere) 

coil 


Figure D. 1 

Comparison of the Calculated and Measured Magnetic Flux Density 
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Appendix E 

Calibration Curve for the Rotameter 


In order to obtain accurate readings, the rotameter used in the experiments was 
calibrated for the test fluid as shown in Figure E. 1 . 



Figure E. 1 

Calibration Curve for the Rotameter 
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